Astronomy & Astrophysics manuscript no. paper 


©ESO 2012 


August 31, 2012 





Dipolar versus multipolar dynamos: the influence of the 
background density stratification 

T. Gastine, L. Duarte, and J. Wicht 

Max Planck Institut fiir Sonnensystemforschung, Max Planck StraBe 2, 37191 Katlenburg Lindau, Germany e-mail: 
gastine@mps . mpg . de 

August 31, 2012 

ABSTRACT 

Context. Dynamo action in giant planets and rapidly rotating stars leads to a broad variety of magnetic field geometries including small 
scale multipolar and large scale dipole-dominated topologies. Previous dynamo models suggest that solutions become multipolar once 
inertia becomes influential. Being tailored for terrestrial planets, most of these models neglected the background density stratification. 
Aims. We investigate the influence of the density stratification on convection-driven dynamo models. 

Methods. Three-dimensional nonlinear simulations of rapidly rotating spherical shells are employed using the anelastic approximation 
to incorporate density stratification. A systematic parametric study for various density stratifications and Rayleigh numbers at two 
different aspect ratios allows to explore the dependence of the magnetic field topology on these parameters. 

Results. Anelastic dynamo models tend to produce a broad range of magnetic field geometries that fall on two distinct branches 
with either strong dipole-dominated or weak multipolar fields. As long as inertia is weak, both branches can coexist but the dipolar 
branch vanishes once inertia becomes influential. The dipolar branch also vanishes for stronger density stratifications. The reason is 
the concentration of the convective columns in a narrow region close to the outer boundary equator, a configuration that favors non- 
axisymmetric solutions. In multipolar solutions, zonal flows can become significant and participate in the toroidal field generation. 
Parker dynamo waves may then play an important role close to onset of dynamo action leading to a cyclic magnetic field behavior 
Conclusions. These results are compatible with the magnetic field of gas planets that are likely generated in their deeper conducting 
envelopes where the density stratification is only mild. Our simulations also suggest that the fact that late M dwarfs have dipolar or 
multipolar magnetic fields can be explained in two ways. They may differ either by the relative influence of inertia or fall into the 
regime where both types of solutions coexist. 
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1. Introduction 

The magnetic fields of planets and rapidly rotating low-mass 
stars are maintained by magnetohydrodynamic dynamos oper- 
ating in their interiors. Scaling laws for convection-driven dy- 
namos successfully predict the magnetic field strengths for both 
types of objects indicating that similar mechanisms are at work 
(Christensen et al. 2009). Recent observations show a broad 
variety of magnetic field geometries, ranging from large-scale 
dipole-dominated topologies (as on Earth, Jupiter and some 
rapidly rotating M dwarfs, see e.g. Donati et al. 2006; Morin 
et al. 2008) to small-scale more complex magnetic structures 
(e.g. Donati et al. 2008; Morin et al. 2010). Explaining the dif- 
ferent field geometries remains a prime goal of dynamo the- 
ory. Unfortunately, the extreme parameters of planetary and 
stellar dynamo regions can not directly be adopted in the nu- 
merical models so that scaling laws are of prime importance 
(Christensen 2010). 

Global simulations of rapidly rotating convection success- 
fully reproduce many properties of planetary dynamos (e.g. 
Christensen & Aubert 2006). They show that the ordering influ- 
ence of the Coriolis force is responsible for producing a dom- 
inant large scale dipolar magnetic field (for a stellar applica- 
tion, see also Brown et al. 2010) provided the Rayleigh num- 
ber is not too large. Multipolar geometries are assumed once 
the Rayleigh number is increased beyond a value where in- 
ertial forces become important (Kutzner & Christensen 2002; 



Sreenivasan & Jones 2006). Christensen & Aubert (2006) sug- 
gest that the ratio of (nonlinear) inertial to Coriolis forces can 
be quantified with what they call the "local Rossby number" 
Ro^ = Urms/^^, where £ is the typical lengthscale of the convec- 
tive flow. Independently of the other system parameters, the tran- 
sition between dipole-dominated and multipolar magnetic fields 
always happens around Ro^ ^0.1. Scaling laws then, for exam- 
ple, predict a dipole-dominated field for Jupiter and a multipolar 
field for Mercury (Olson & Christensen 2006). The latter may 
explain the strong quadrupolar component in the planet's mag- 
netic field (Christensen 2006). Earth lies at the boundary where 
the field is dipolar most of the time but occasionally forays into 
the multipolar regime allowing for magnetic field reversals. 

The scaling laws based on Boussinesq simulations are geared 
to model the dynamo in Earth's liquid metal core where the 
background density and temperature variations can be neglected. 
Their application to gas giants and stars where both quantities 
vary by orders of magnitude is therefore questionable (Chabrier 
& Kiiker 2006; Nettelmann et al. 2012). Recent compressible 
dynamo models of fully-convective stars suggest a strong influ- 
ence of the density stratification on the geometry of the mag- 
netic field: while the fully compressible models of Dobler et al. 
(2006) have a significant dipolar component, the strongly strati- 
fied anelastic models of Browning (2008) tend to produce mul- 
tipolar dynamos (see also Bessolaz & Brun 2011). These diflTer- 
ences stress the need of more systematic parameter studies to 
clarify the influence of density stratification. 
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In addition, most of the previous Boussinesq studies have 
employed rigid flow boundary conditions for modelling terres- 
trial dynamo models. Stress-free boundary conditions, more ap- 
propriate for gas planets and stars, show a richer dynamical be- 
haviour, including hemispherical dynamos and bistability where 
dipole-dominated and multipolar dynamos coexist at the same 
parameters (e.g. Grote & Busse 2000; Busse & Simitev 2006; 
Goudard & Dormy 2008; Simitev & Busse 2009; Sasaki et al. 
2011). This challenges the usefulness of the Ro^ criterion since 
bistable cases exist significantly below Ro^ ^0.1 (e.g. Schrinner 
et al. 2012). 

Here we adopt the anelastic approximation (e.g. Gilman & 
Glatzmaier 1981; Braginsky & Roberts 1995; Lantz & Fan 1999) 
to explore the eflTect of a background density stratification on the 
dynamo process while filtering out fast acoustic waves. We con- 
duct an extensive parameter study where we vary the degree of 
stratification and the Rayleigh number to determine the param- 
eter range where dipole-dominated fields can be expected. The 
anelastic approximation and the numerical methods are intro- 
duced in section 2. In section 3, we present the results of the 
parametric study and describe the diflTerent dynamo regimes be- 
fore relating our results to observations in the concluding sec- 
tion 4. 



2. The dynamo model 

2.1. An anelastic formulation 

We consider MHD simulations of a conducting ideal gas in a 
spherical shell rotating at a constant frequency Q about the z 
axis. Convective motions are driven by a fixed entropy contrast 

between the inner radius and the outer radius r^. Following 
the previous parametric studies of Christensen & Aubert (2006), 
we use a dimensionless formulation, where the shell thickness 
d = Vg- Vi is the reference lengthscale, the viscous diflTusion time 
(f Iv is the reference timescale and ^pto^jiAD. is the magnetic 
scale. Density and temperature are non-dimensionalised using 
their values at the outer boundary ptop and Ttop and A^ serves 
as the entropy scale. The kinematic viscosity v, magnetic diflTu- 
sivity A, magnetic permeability yu, thermal diflfusivity k and heat 
capacity Cp are assumed to be constant. 

Following the anelastic formulation by Gilman & Glatzmaier 
(1981), Braginsky & Roberts (1995) and Lantz & Fan (1999) we 
adopt a non-magnetic, hydrostatic and adiabatic background ref- 
erence state that we denote with overbars in the following. It is 
defined by the temperature profile dT/dr = -g/Cp and the as- 
sumption of a poly tropic gas p(r) = T"^, where m is the poly- 
tropic index. Gravity typically increases linearly with radius in 
Boussinesq models where the density is approximately constant. 
Many anelastic models on the other hand assume that the den- 
sity is predominantly concentrated below the dynamo zone so 
that g oc 1/r^ provides a better approximation (e.g. Glatzmaier 
1984; Jones et al. 2011). We adopt the form 
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Here, r] = Vilvo is the aspect ratio of the spherical shell and Np 
is the number of density scale heights covered by the density 
background. 

Under the anelastic approximation the dimensionless equa- 
tions governing convective motions and magnetic field genera- 
tion are then given by 
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where u, B,p and s are velocity, magnetic field, pressure and en- 
tropy disturbance of the background state, respectively. S is the 
traceless rate-of- strain tensor with a constant kinematic viscosity 
given by 
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where dij is the identity matrix. In addition to the aspect ra- 
tio of the spherical shell r] and the two parameters involved in 
the description of the reference state {Np and m), the system of 
equations (4-7) is governed by four dimensionless parameters, 
namely the Ekman number E, the Prandtl number Pr, the mag- 
netic Prandtl number Pm and the Rayleigh number Ra defined 
by 



Q.d^ K A CpVK 

where gtop is the reference gravity at the outer boundary r^. 
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with constants gi and g2 to allow for both alternatives which then 
leads to the following background temperature profile (see also 
Jones et al. 2011; Gastine & Wicht 2012) 
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2.2. The numerical method 

The numerical simulations in this parameter study have been 
computed with the anelastic version of the code MagIC (Wicht 
2002; Gastine & Wicht 2012), which has been vaHdated in 
an anelastic dynamo benchmark (Jones et al. 2011). To solve 
the system of equations (4-7), pu and B are decomposed into 
poloidal and toroidal contributions 



pu = Vx(VxWer)-\-VxZer , B = Vx(VxC^^)-fVxA^^. (10) 
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The unknowns ly, Z, A, C, ^ and p are then expanded in spher- 
ical harmonic functions up to degree ^max in colatitude 6 and 
longitude and in Chebyshev polynomials up to degree A^^ in 
the radial direction. An exhaustive description of the complete 
numerical method and the associated spectral transforms can be 
found in (Gilman & Glatzmaier 1981). Typical numerical resolu- 
tions employed in this study range from (A^^ = 61, ^max = 85) for 
Boussinesq simulations close to onset to (A^^ = 81, ^max = 133) 
for the more demanding cases with a significant density contrast. 

In all the simulations presented in this study, we have as- 
sumed constant entropy boundary conditions at and r^. The 
mechanical boundary conditions are either no slip at the inner 
and stress-free at the outer boundary or stress-free at both limits 
(see Tab. 2). The magnetic field is match to a potential field at 
both boundaries. 



2.3. Non-dimensional diagnostic parameters 

To quantify the impact of the diff'erent control parameters on 
magnetic field and flow, we analyse several diagnostic non- 
dimensional properties. The typical rms flow amplitude in the 
shell is either given as the magnetic Reynolds number Rm = 
Uvm^djA or the Rossby number Ro = u^^^jVld. Rm is a mea- 
sure for the ratio of magnetic field generation to Ohmic dis- 
sipation and thus an important quantity for any dynamo. The 
Rossby number quantifies the ratio between inertia and Coriolis 
forces but Christensen & Aubert (2006) demonstrated that the 
local Rossby number Ro^ = Wrms/^^ is a more appropriate mea- 
sure, at least concerning the impact of inertia on the magnetic 
field geometry. The typical flow lengthscale I can be calculated 
based on the spherical harmonic flow contributions 



with 
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Here, is the velocity field at a given spherical harmonic degree 
i and (• • • ) corresponds to an average over time and radius. 

The magnetic field strength can, for example, be measured 
by the Elsasser number A = B^^/pjio^^, which is supposed to 
provide an estimate for the ratio of Lorentz to Coriolis forces. 
This estimate may be rather far from the true force balance (e.g. 
Wicht & Christensen 2010). The modified Elsasser number sug- 
gested by Soderlund et al. (2012) and defined by 
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provides a more appropriate measure, being more directly re- 
lated to the ratio of the two respective terms in the Navier-Stokes 
equation (4) when approximating the lengthscale entering in the 
curl in the Lorentz force by In addition to the typical spher- 
ical harmonic degree ^5, the typical order is also used since 
both seem to contribute (Soderlund et al. 2012): 
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Table 1. Critical Rayleigh numbers and corresponding az- 
imuthal wave numbers for the simulations considered in this 
study (for all of them E = lO'^ and Pr = 1). 
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Finally, the geometry of the magnetic field is quantified by 
its dipolarity 



/dip — 



) 5^^^ sin 6 do d(p 
§B^smOded(f) 



(15) 



that measures the ratio of the magnetic energy of the dipole to 
the total magnetic energy at the outer boundary r^. 

2.4. A parameter study 

In all the simulations presented in this study, the Ekman num- 
ber is kept fixed to a moderate value of E = 10""^, which al- 
lows to study a large number of cases. The Prandtl number is 
set to 1 and the magnetic Prandtl number to 2. Following Jones 
& Kuzanyan (2009) and our previous hydrodynamical models 
(Gastine & Wicht 2012), we adopt a polytropic index m = 2 
for the reference state. We consider two diflferent aspect ratios: 
thick shells with r] = 0.2 and thin shells with r] = 0.6. Dynamo 
models in thin shells rely on the same setup as our previous hy- 
drodynamical study (Gastine & Wicht 2012). Gravity is adapted 
to better reflect the difl'erent shells, i.e. we use linear gravity in 
the thick shells (gi = I, g2 = 0) and g oc 1/r^ in the thin shells 
(gi = 0, g2 = 1). Both cases also diff'er in the flow boundary 
conditions: thin shell cases assume stress-free conditions at both 
boundaries, thick shell cases assume a rigid inner boundary. 

Being mostly interested in the eff'ects of the density strat- 
ification, we have performed numerical simulations with dif- 
ferent density contrasts spanning the range from Np = (i.e. 
Boussinesq) to Np = 3 (i.e. Pbot/Ptop - 20). For each den- 
sity stratification, we vary the Rayleigh number from simula- 
tions close to onset of dynamo action to 10-50 times the criti- 
cal Rayleigh number Ra^ for the onset of convection. Ra^ varies 
with Np and we list respective values for the cases explored here 
in Tab. 1 . These have been determined with the anelastic linear 
stability analysis code developed by Jones et al. (2009). 

As will be shown below, we find several cases of bistability 
where a strong dipole-dominated and a weak multipolar solu- 
tion coexist at identical parameters. The starting solution for our 
time integration then determines which of the two solutions the 
simulation will assume. Most of the numerical simulations there- 
fore have been initiated with both a strong dipolar magnetic field 
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Fig. 1. Relative dipole strength plotted against the local Rossby 
number Ro^ for two different aspect ratios: rj = 0.2 (upper panel) 
and 77 = 0.6 (lower panel). Each type of symbol is associated to 
a given density stratification (from Boussinesq, i.e. A/p = to 
stratified simulations corresponding to Np = 3). The size of the 
symbol has been chosen according to the value of the modified 
Elsasser number (Eq. 13). Vertical lines are tentative transitions 
between dipolar and multipolar dynamos. 



(with A ~ 1) and a weak multipolar magnetic field (model names 
end with a 'd' or a 'm' accordingly in Tab. 2). Altogether, more 
than 100 simulations have been computed, each running several 
magnetic diffusion time to ensure that a stable magnetic config- 
uration has been reached. Table 2 lists the input parameters of 
all cases along with key solution properties. 

3. Results 

3.1. Dynamo regimes 

Figure 1 shows how the dipolarity /dip depends on the local 
Rossby number in the thick shell (top) and the thin shell (bottom) 
simulations. In both geometries we find two distinct branches: 
The upper branch corresponds to the dipolar regime at /dip > 0.5 
with strong magnetic fields at modified Elsasser number be- 
tween = 0.1 and A^ = 1. This indicates a strong Lorentz 
force so that these dynamos likely operate in the so-called mag- 
netostrophic regime where Coriolis forces, pressure gradients. 
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Fig. 2. Regime diagram for dynamo simulations as a function of 
supercriticality and density stratification for two diflTerent aspect 
ratios: t] = 0.2 (upper panel) and t] = 0.6 (lower panel). Red cir- 
cles correspond to dipolar dynamos, blue squares to multipolar 
dynamos and black crosses are decaying dynamos. The size of 
the symbols depend on the value of the modified Elsasser num- 
ber (Eq. 13). The grey-shaded area highlights the dipolar region, 
while the dashed lines correspond to the critical Ro^c values that 
mark the limit between dipolar and multipolar dynamos (see the 
vertical lines in Fig. 1). 



buoyancy and Lorentz forces contribute to the first order Navier- 
Stokes equation. In contrast, the models belonging to the lower 
branch have a multipolar field geometry at /dip < 0.2 and much 
weaker magnetic field with A^ < 5 x 10"^. Since the Lorentz 
force may thus not enter the first order Navier- Stokes equation 
these cases are likely geostrophic. 

For both aspect ratios, the dipolar branch is bounded by a 
critical local Rossby number Ro^^ beyond which all cases are 
multipolar. In the thin shell with 7] = 0.6, the value Ro^^ - 0.15 
is slightly larger than the value of Ro^c - 0-12 predicted by 
Christensen & Aubert (2006) in their geodynamo models with 
J] = 0.35. In the thick shell with 7] = 0.2, it is smaller at 
Ro£c - 0.08. The critical local Rossby number thus seems to in- 
crease with the inner core size. A similar dependence on the as- 
pect ratio has been already found by Aubert et al. (2009). In con- 
trast to the previous Boussinesq studies that used rigid boundary 
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conditions (e.g. Christensen & Aubert 2006; Aubert et al. 2009), 
the multipolar branch also extends here below the critical Ro^c 
where the dipolar and the multipolar branch now coexist. This 
multipolar dynamo branch for Ro^ < 0.1 has also been found 
in the Boussinesq models of Schrinner et al. (2012) with stress- 
free boundary conditions. We further discuss the bistability phe- 
nomenon in section 3.2. 

Figure 2 illustrates how the field morphology depends on 
the stratification and on the supercriticality Ra/Ra^. The value 
Np ^ 1.8 marks an important regime boundary: below this 
value most cases are bistable (nested circles and squares), while 
above this value only multipolar cases remain (squares only). 
The reason for this boundary will be considered in section 3.3. 
The dashed black lines in Fig. 2 mark the Ro^c boundary which 
moves to lower supercriticalities for increasing density stratifi- 
cations. Larger stratifications promote larger flow amplitudes as 
well as smaller lengthscales since convection is progressively 
concentrated in a narrowing region close to the outer boundary 
(see Table 1, Fig. 9 and Jones & Kuzanyan 2009; Gastine & 
Wicht 2012). Both eflfects lead to larger Ro^ for a given super- 
criticality. 

For Np < 1.8 the boundary for onset of dynamo action 
also seems to move to smaller supercriticalities when Np in- 
creases. The reason is a mixture of a decreasing critical mag- 
netic Reynolds number and an increasing flow amplitude for a 
given supercriticality. Beyond Np = 1.8, the critical magnetic 
Reynolds number once more increases and the onset of dynamo 
action moves to larger supercriticalities. The onset of dynamo 
action and the Ro^^ boundary are less aflfected by the stratifi- 
cation for the thinner shell. One reason may be that the flow 
lengthscale is already quite small for thinner shells even in the 
Boussinesq case (see Tab. 1 and Al-Shamali et al. 2004). The dif- 
ferences in the gravity profile and the flow boundary condition 
may also play a role here. 

The dipolar branch reaches down to lower local Rossby num- 
bers than the multipolar branch (see Fig. 1). This is where we 
find examples of subcritical dynamo action (models Id, 6d and 
12d in Tab. 2), where the dynamo is only successful when started 
with a sizable magnetic field strength (Morin & Dormy 2009). 
Simulations started with a weak field do not have the option to 
run to the multipolar branch and the field simply decays away. 

3.2. Bistability 

As illustrated on Fig. 2, most of our models at Np < 1.8 show 
bistability. Exceptions are the subcritical thick shell cases and 
models 13d/m and 14d/m (see Tab. 2), where even small ini- 
tial magnetic fields developed to strong dipoles. The top panel 
of Fig. 3 shows the dipole-dominated solution with, however, 
stronger dynamo action in the southern hemisphere. The lower 
panel of Fig. 3 illustrates the weak multipolar solution, which 
has pronounced equatorially symmetric and non-axisymmetric 
components (this is a feature also visible in the more stratified 
cases shown in the following Figs. 7-8). The dipolar solution 
has a three times larger modified Elsasser number. Generally, 
is twice to ten times larger in the dipole-dominated than 
in the multipolar counterpart of bistable cases (see Fig. 1 and 
Tab. 2 for further details). This suggests that the two dynamo 
branches are also characterised by diflTerent force balances, at 
least in the bistability region (i.e. Ro^ < 0.1): the dipolar branch 
is magnetos trophic, while the multipolar one is more likely to 
be geostrophic. Beyond Ro^ > 0.1, however, the multipolar dy- 
namos can be strong enough to yield a magneto strophic force 
balance due to the larger Rayleigh numbers (Brun et al. 2005). 




Fig. 3. Left panels: iso-contours of axisymmetric azimuthal field 
and poloidal field lines. Right panels: surface radial field Br(r = 
To). From a model with 77 = 0.2, A/p = 1.7, Ra = 9 x 10^ initiated 
with a dipolar magnetic field (upper panels) and a multipolar 
one (lower panels). Red (blue) correspond to positive (negative) 
values. 
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Fig. 4. Time-averaged kinetic and magnetic energy spectra plot- 
ted against the spherical harmonic degree i. From a model with 
7] = 0.2, A/p = 1.7, Ra = 9 X 10^ initiated with a dipolar magnetic 
field (black lines) and a multipolar one (blue lines), correspond- 
ing to Fig. 3. 



Figure 4 shows the corresponding time-averaged kinetic and 
magnetic spectra for the two bistable solutions displayed in 
Fig. 3. The two kinetic energy spectra nearly coincide, emphasis- 
ing the similarities of the flow in the two solutions. We observe 
a broad plateau for degrees 4 < ^ < 30 that reflects the convec- 
tive driving (for this model, the critical azimuthal wavenumber 
at onset of convection is m = 6, see Tab. 1). The magnetic spec- 



5 



T. Gastine et al.: Dipolar versus multipolar dynamos: the influence of the background density stratification 



0.1 



o 



0.01 



10- 



1) Strong dipole 
Weak shear 




^Np = 0.5 

^p=i 



□ ^^P = 2 
^ Np = 2.5 




2) Weak dipole 
Strong shear 



3) Weak dipole 
Weak shear 



-0, Ra= 1.0 X 10^ 



-3, Ra= 1.1 X 10^ 



0.01 



0.03 



0.1 

Etor axi/Etor 



0.3 



Fig. 5. Relative dipole strength plotted against the ratio of ax- 
isymmetric toroidal kinetic energy divided by toroidal energy. 
Red (grey) symbols correspond to simulations in thick (thin) 
shells (77 = 0.2 and 77 = 0.6, respectively). Each type of symbol is 
associated with a given density stratification (from Boussinesq, 
i.e. Np = to Np = 3). The sizes of the symbols have been 
chosen according to the values of the modified Elsasser number 
(Eq. 13). 



tra show the diflTerences already visible in the surface magnetic 
fields displayed in Fig. 3. The amplitude of the dipolar compo- 
nent is more than one order of magnitude lower, higher-order 
components drop by roughly 50%. A mild maximum for ^ = 4 
corresponds to the large-scale contribution visible in Fig. 3. For 
spherical harmonic degrees ^ > 30, the spectra follow some clear 
power-law behaviour. Both magnetic and kinetic energy spectra 
have a steep decrease (between and very similar to pre- 
vious stellar dynamo models (e.g. Brun et al. 2005; Browning 
2008) or quasi-geostrophic kinematic dynamos (e.g. SchaeflTer 
& Cardin 2006). 

The coexistence of a dipolar and a multipolar branch has 
already been reported for Boussinesq models with stress-free 
(Busse & Simitev 2006; Schrinner et al. 2012) or mixed mechan- 
ical boundary conditions (Sasaki et al. 201 1). For rigid boundary 
conditions only one such case has been found by Christensen 
& Aubert (2006). Stress-free boundaries allow stronger zonal 
winds to develop which seem to play a key role here. These 
geostrophic flows (constant on coaxial cylinders) are maintained 
by Reynolds stresses, i.e. by the correlation between the zonal 
and the cylindrically radial Us velocity components (e.g. Gastine 
& Wicht 2012). For stress-free boundary conditions their ampli- 
tude is limited by the weak bulk viscosity. When rigid boundary 
conditions are employed, the much stronger boundary friction 
severely brakes these zonal winds. 

Figure 5 displays the dipolarity of the surface field against 
the ratio of axisymmetric toroidal to total toroidal kinetic en- 
ergy, which is a good proxy of the relative amplitude of zonal 
winds. Strongly dipolar solutions cluster in the upper left corner 
where zonal winds are weak. For strong zonal winds only weak 
multipolar solutions can be found. For bistable cases, one solu- 
tion belongs to the first type while the other one belongs to the 
second category. There clearly is a competition between strong 
zonal winds and strong dipolar magnetic fields. The third type 
of multipolar and sometimes strong solutions with mostly weak 
zonal winds are mainly strongly stratified cases that we will dis- 




Q-effect 



Fig. 6. Role of the D-eflfect in the production of the toroidal field 
for two selected simulations in thick shells (upper panels) and 
in thin shells (lower panel). For each simulation, the Il-eflfect 
(color levels) and the axisymmetric azimuthal field (solid and 
dashed lines) are displayed on the left part, while the axisymmet- 
ric zonal flow is given on the right part. Red (blue) correspond 
to positive (negative) values. 



cuss in section 3.3. The strong multipolar fields can have a strong 
impact on the zonal flows, in a similar way to the stellar models 
of Brun et al. (2005). 

The stronger zonal winds in the multipolar cases go along 
with a change in dynamo mechanism which is illustrated in 
Fig. 6. In the dipolar cases (left panels), zonal flows are weak 
and mainly driven by thermal wind eflTects. The stronger Lorentz 
force can balance the Coriolis force and allows for these strongly 
non-geostrophic motions (Aubert 2005). The Q-eflTect, i.e. the 
production of axisymmetric toroidal magnetic field by zonal 
wind shear, plays only a secondary role (left half of left pan- 
els in Fig. 6) so that the dynamo is of the type in the mean 
field nomenclature (Olson et al. 1999), similar to the planetary 
dynamo models by Christensen & Aubert (2006) or mean-field 
models of fully convective stars (Chabrier & Kiiker 2006). In the 
multipolar cases (right panels in Fig. 6), significant zonal winds 
driven by Reynolds stresses develop close to the outer boundary 
and the associated D-eflTect plays an important role for toroidal 
magnetic field production. Dynamos on the multipolar branch 
(for Np < 2) are thus of the afl or a^D. type. Many of these 
solutions show distinct oscillations connected to Parker dynamo 
waves that we discuss further in section 3.4. 

3.3. Influence of stratification 

For Np > 1.8 only the multipolar branch remains. Figures 7 and 8 
illustrate how the convective flow and the surface magnetic field 
evolve on increasing the density stratification for simulations 
with very similar local Rossby numbers. All these models have 
been initiated with a strong dipolar field. For both aspect ratios. 
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Np = 0, Ra = 2.0 X 10^ Np = 1, Ra = 1.0 x 10^ Np = 2, Ra = 1.0 x 10^ Np = 3, Ra = 1.1 x 10^ 




Fig. 7. Surface radial field Br(r = ro) (upper panels) and radial velocity close to the outer boundary Ur{r = 0.9 r^,) (lower panels) for 
four numerical simulations in thick shells (i.e. t] = 0.2) with similar local Rossby number Ro^ ^ 0.04 - 0.05 (see Tab. 2). Red (blue) 
correspond to positive (negative) values. 



Np = 10-2, Ra = 2.0 X 10^ Np = 1, Ra = 1.5 x 10^ Np = 2, Ra = 2.0 x 10^ Np = 3, Ra = 3.0 x 10* 




Fig. 8. Surface radial field Br(r = ro) (upper panels) and radial velocity close to the outer boundary Ur(r = 0.9 r^,) (lower panels) for 
four numerical simulations in thin shells (i.e. t] = 0.6) with similar local Rossby number Ro^ ^ 0.07 - 0.08 (see Tab. 2). Red (blue) 
correspond to positive (negative) values. 



the weakly stratified cases (i.e. Boussinesq and Np = I) have 
a dipole-dominated magnetic field, while the stronger stratified 
cases (i.e. Np = 2 and Np = 3) have multipolar fields with signif- 
icant non-axisymmetric contributions. As a consequence of the 
background density contrast, the typical lengthscale of convec- 
tion gradually decreases with Np (see also Tab. 1 and Gastine 
& Wicht 2012). While the imprints of these smaller convective 
scale are visible in the surface magnetic fields, the dynamos are 
always dominated by larger scale features even in the strongly 
stratified cases. For example, for the two models at Np = 3, a 
clear wave number m = 2 signature emerges similar to the pat- 
terns already observed in the Boussinesq models of Goudard & 
Dormy (2008). This magnetic mode can only be clearly identi- 
fied for moderate magnetic Reynolds numbers Rm < 200 close 
to the onset of dynamo action. For larger Rm, smaller scale more 
complex features once more take over. 

Some models also show a stronger concentration of the mag- 
netic field in one hemisphere. Similar eflfects can be found for 
some multipolar simulations at low stratifications (see section 
3.4). Grote & Busse (2000) and Busse & Simitev (2006) report 
that hemispherical dynamo action is typical for simulations with 
stress-free boundaries and Prandtl numbers around unity. 

The collapse of the dipolar branch for A/p > 2 is caused 
by the strong concentration of convective features close to 



the outer boundary where density decreases most drastically. 
Figure 9 illustrates the evolution of the convective columns when 
Np increases while the local Rossby numbers remains similar 
(see Figs. 7-8 for the corresponding magnetic fields and ra- 
dial flow structures). As it has already been demonstrated by 
non-magnetic simulations (Jones et al. 2009; Gastine & Wicht 
2012), convection first develops close to the inner boundary for 
Boussinesq and weakly stratified models (Np < 1). When in- 
creasing Np, however, the convective columns gradually move 
outward and become confined to a thin region close to the equa- 
tor (Np > 2). This goes along with a decrease of the typical 
lengthscale (see also Tab. 1 for the critical azimuthal wave num- 
bers at onset). According to Gastine & Wicht (2012), this is a 
consequence of the fact that buoyancy and thus the eflTective lo- 
cal Rayleigh number becomes much larger at the outer boundary 
than in the interior for strongly stratified models. 

The kinetic helicity *K = w'-Vxw'isa key ingredient 
in the induction process via the Of-eflfect (e.g. Moflfatt 1978). 
Figure 10 illustrates how the helicity increasingly concentrates 
closer to the outer boundary equator when the stratification in- 
tensifies simply because the convective columns are the main 
carriers of helicity. Mean field dynamo models have demon- 
strated that a similar concentration of the Of-eflfect can promote 
non-axisymmetric dynamo modes of low spherical harmonic or- 
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iVp = 0, Ra = 2 X 10^ 



= 1, Ra = 10^ 



A^. = 2, Ra = 10^ 




Np = 10-2, Ra = 2 X 10^ 




Np = 1, Ra = 1.5 X 10' 




Np = 2, Ra = 2 X 10^ 




Np = 3, Ra= 1.1 X 10^ 




No = 3, Ra = 3 X 10^ 




Fig. 9. Isocontours and equatorial cut of vorticity along the axis of rotation cl)^ = (V x m)^ for different density stratification in thick 
shells (7/ = 0.2, upper panels) and in thin shells (7/ = 0.6, lower panels). These models correspond to the cases displayed in Figs. 7-8. 
Positive (negative) values are rendered in red (blue). For each model, the white sphere corresponds to the inner radius. 



Np = 




t = O.OOe+OO t = 2.25e-02 t = 4.72e-02 t = 7.22e-02 t = 9.65e-02 







t = O.OOe+00 t = 5.97e-02 t = 1.20e-01 t = 1.80e-01 t = 2.40e-01 




Fig. 10. Azimuthal average of kinetic helicity • V x for the 
models displayed in Fig. 9. Positive (negative) values are ren- 
dered in red (blue). 



der (typically m ~ 1 - 2, see Riidiger et al. 2003; Jiang & Wang 
2006). Non-axisymmetric a^D. mean-field models also suggest 
a lower critical dynamo number for non-axisymmetric dynamo 
modes (e.g. Ruzmaikin et al. 1988; Bassom et al. 2005). These 
findings may explain the onset of the new larger wave number 
magnetic mode in our strongly stratified cases. 



Fig. 11. Time evolution of axisymmetric toroidal magnetic field 
B(p for a simulation in a thick shell (t] = 0.2, Np = 3, Ra = 
1.1 X 10^, upper panel) and a simulation in a thin shell (t] = 0.6, 
A/p = 1, Ra = 10^, lower panel). Times are given here in mag- 
netic diffusion units. Red (blue) correspond to positive (negative) 
values. 



3.4. Parker waves 

The aQ. or a^Q dynamos on the multipolar branch also show a 
distinct and interesting oscillatory time dependence. Figures 11- 
12 illustrate two examples. The magnetic field is weak, multi- 
polar and dominated by non-axisymmetric components in both 
cases (see for example Figs. 7-8). For the thin shell (lower pan- 
els), a dipolar solution coexists at identical parameters. Toroidal 
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90° 






r = O.y 



5.0 5.5 

Time 



b) 



-90° 



r = 0.9 



Time 



Fig. 12. Butterfly diagrams: time evolution of axisymmetric 
toroidal magnetic field 5^ at r = 0.9ro for the two models of 
Fig. 11. Time is given here in magnetic diflTusion units. Red 
(blue) correspond to positive (negative) values. 

magnetic field first emerges at low latitudes, then travels towards 
the poles and finally vanishes roughly where the tangent cylinder 
seems to prevent a further movement towards even higher lati- 
tudes (Schrinner et al. 2011). When the patches have travelled 
half way, the cycle starts over with field of opposite polarity ap- 
pearing close to the equatorial plane. The thick shell solution 
(top panels) demonstrates that during some periods one hemi- 
sphere clearly dominates. The characteristic period is roughly 
O.lr^ - 0.2t^ where = d^/A is the magnetic diflTusion time. 
These are Parker dynamo waves, as we will demonstrate below, 
very similar to those that have been previously identified in some 
Boussinesq simulations (Goudard & Dormy 2008; Schrinner 
et al. 2011; Simitev & Busse 2012). 

These coherent oscillations lead to the periodic patterns ob- 
served in the "butterfly-diagram" style illustrations shown in 
Fig. 12. The diagrams highlight the direction of travel which in 
Parker waves is controlled by the gradient of zonal flows: they 
travel polewards when like in our simulations the diflTerential ro- 
tation decreases with depth (e.g. Yoshimura 1975). 

The mean-field formalism allows to derive a dispersion rela- 
tion for these dynamo waves (Busse & Simitev 2006; Schrinner 
et al. 2011). In this formalism, magnetic field and flow are de- 
composed into the axisymmetric or mean contributions denoted 
by overbars 



B = Bp + Beef, = V X (A ^0) Beef, ; u = u( 



(16) 



and the primed non-axisymmetric or fluctuating contributions B^ 
and u\ This allows to reformulate the axisymmetric part of the 
induction equation (5) to 



— =aB-\- — V^A, 
dt Pm 

dB 1 . 

— = Bp-Vu + (Vx aBp) • + 
\ at y T 



(17) 



where the fluctuating mean electromotive force is approximated 

by an a-eflTect, i.e. by x B^ = aB. 

To simplify this system we assume a homogeneous a and 
also adopt the plane layer formalism introduced by Parker 
(1955), where the cartesian coordinates (x,y,z) correspond to 
(0, 6, r). Further assuming that u depends only on z (radius) leads 
to 



— =aB+ — V^A, 
dt Pm 

dB „9 du OA 1 „9 
dt dz oy Pm 



(18) 



The ansatz A,Boc Qxp(i(kyy -\- k^z) + At) with A = r + ico allows 
to derive the following dispersion relation 



Pm / dz 



(19) 



For aQ. dynamos we concentrated on in the following, the imag- 
inary part of A provides 



1 du ^ 



1/2 



(20) 



To further simplify this dispersion relation we assume that z vari- 
ations are of the order of the shell radius, i.e. ky ~ l/r^ and that 
the shear can be approximated by du/dz ~ u^/d ~ RCzon- 

The a is particularly difl&cult to estimate. A full derivation 
of the a tensor would, for example, require to employ the test- 
field method (e.g. Schrinner et al. 2007). For homogeneous and 
isotropic MHD turbulence, however, it may at least crudely be 
approximated via the fluctuating kinetic helicity: 



(21) 



where is the typical lifetime of a convective feature (e.g. 
Brandenburg & Subramanian 2005). Following Brown et al. 
(2010), we use = Hp/u' when the local density scale height 
Hp = (dlnp/dr)~^ is smaller than d and = d/u' otherwise. 
This finally leads to the simplified dispersion relation for Parker 
waves 



<jj 



oRezo 

2 To 



1/2 



(22) 



Figure 13 shows the comparison between the frequencies 
predicted by this simplified dispersion relation and those found 
in the numerical models. The later frequencies are obtained us- 
ing a Fourier transform of the butterfly diagrams displayed in 
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30 



10 
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0^p = 3 



10 30 100 



300 



Fig. 13. Frequencies of oscillatory dynamos plotted against the 
dispersion relation of Parker waves given in Eq. (22). Red (grey) 
symbols correspond to simulations in thick (thin) shells (t] = 0.2 
and 7] = 0.6, respectively). Each type of symbol is associated to 
a given density stratification. 



Fig 12. The resulting B^piO^oj) is then integrated over the co- 
latitude 6 to derive the power spectrum. Even though individ- 
ual predictions fail by up to a factor two, the general agreement 
is quite convincing given the numerous approximation involved 
and is similar to that reached in previous Boussinesq studies 
(Busse & Simitev 2006; Schrinner et al. 2011). This confirms 
that the oscillations in our multipolar cases are indeed Parker 
waves. When the magnetic Reynolds number becomes too large 
(Rm > 150 - 200), the coherence of these oscillations is gradu- 
ally lost and a definite frequency is increasingly difliicult to de- 
termine. This limits our analysis to lower RCzon and a values and 
to frequencies cl) < 300. 

4. Discussion and conclusions 

We have investigated the influence of background density strat- 
ification on convection-driven dynamos in a rotating spherical 
shell. The use of the anelastic approximation allowed us to ex- 
clude sound waves and the related short time steps (Gilman & 
Glatzmaier 1981; Clune et al. 1999; Jones & Kuzanyan 2009). 
Previous Boussinesq results have shown that inertial eflTects play 
a decisive role for determining the magnetic field geometry. 
When inertia becomes influential, only multipolar solutions with 
weak magnetic fields are possible. When inertia is weak, two 
solutions can coexist: dipole-dominated solutions with strong 
magnetic fields are found for stress-free as well as rigid bound- 
ary conditions. For stress-free and mixed boundary conditions, a 
multipolar branch is found at identical parameters (e.g. Simitev 
& Busse 2009; Schrinner et al. 2012). Alternatively, the recent 
study by Soderlund et al. (2012) suggests that the transition be- 
tween dipolar and multipolar dynamos may occur when inertia 
becomes larger than viscous forces. This transition is accompa- 
nied by an abrupt decrease of the kinetic helicity. These results 
however seem to be in contradiction with previous studies where 
inertia is always larger than viscosity (e.g. Wicht & Christensen 
2010) and further investigations are required to clarify this con- 
tradiction. 

Our anelastic simulations confirm this scenario for mild 
stratifications corresponding to A^p < 1.8. The reason for the 



bistability is a competition between zonal winds and dipolar 
magnetic fields. Strong dipolar magnetic fields prevent signifi- 
cant zonal winds to develop. Strong zonal winds, on the other 
hand prevent the production of significant dipolar fields. The 
two branches also diflTer in the induction mechanism. Strong 
zonal winds promote an D-eflTect which leads to a aQ. or a^D. 
type of dynamo, while dipole-dominated magnetic fields are 
typically generated in an process. The sizable axisymmet- 
ric toroidal fields produced by an D-eflfect typically lead to a 
coherent cyclic time evolution of the magnetic field for moder- 
ate magnetic Reynolds numbers (Rm < 200). This is consistent 
with previous Boussinesq studies (e.g. Schrinner et al. 2007) and 
numerical models of young solar-type stars (Brown et al. 2011) 
and has been identified as Parker waves. Contrary to what is ob- 
served in the solar cycle, these waves start at the equator and 
travel towards the poles because of the opposite sign in the zonal 
shear. 

For stronger stratification with Np > 1.8, the dipolar branch 
is lost and close to onset of dynamo action a new magnetic mode 
characterised by a large wave number m = 2 appears. This is 
likely due to the concentration of the a-eflTect into a narrow 
region close the equator. According to mean-field models, this 
would preferentially promote non-axisymmetric dynamos (see 
also Chabrier & Kiiker 2006) with large wave numbers. The 
collapse of the dipolar branch may explain the diflTerences be- 
tween the weakly stratified and significantly dipolar simulations 
by Dobler et al. (2006) (pbot/Ptop - 5, i.e. Np = 1.6) and the 
strongly stratified anelastic and multipolar models by Browning 
(2008) (Pbot/Ptop - 100, i.e. Np = 4.6). 

Numerical limitations force us to use excessively large dif- 
fusivities in our simulations. Ekman numbers are thus orders of 
magnitude too large and Reynolds number orders of magnitude 
too low. When extrapolated to fast rotating planets and stars this 
type of models nevertheless provides, for example, realistic mag- 
netic field strengths (Olson & Christensen 2006; Christensen 
et al. 2009). This gives us confidence to compare the observed 
field geometries of planets and stars with predictions based on 
our simulation results. 

Our result are compatible with the dipole-dominated mag- 
netic fields on Jupiter and Saturn that are generated in their 
deeper metallic envelopes where the density stratification is only 
mild (roughly A/p ~ 1.5 - 2, see Heimpel & Gomez Perez 2011; 
Nettelmann et al. 2012). Since the local Rossby number is very 
small for these planets, however, bistability seems to be an op- 
tion. This is also the case for Uranus and Neptune, where Ro^ 
and the density contrast within the dynamo region are small 
(Hubbard et al. 1991; Olson & Christensen 2006). Their multi- 
polar magnetic field would then suggest that these planets would 
then occupy the alternative branch offered by the bistability phe- 
nomenon. 

Concerning rapidly rotating low-mass stars that may also fall 
into the low Ro^ regime, the spectropolarimetric observations of 
Morin et al. (2010) suggest that late M stars with very similar pa- 
rameters (mass and rotation rate) come in two categories: some 
stars present a strong dipole-dominated magnetic field while oth- 
ers show weaker and multipolar magnetic structures. These two 
geometries may represent the two coexisting dynamo branches 
at smaller local Rossby numbers. 

For the bistability to be a viable explanation, however, our 
simulations suggest that the dynamos must operate in a region 
with moderate density stratification (Np ~ 1 - 2) and super- 
criticality (Ra/Ra^ ~ 10 - 30). Since the stellar dynamos op- 
erate presumably far from onset of convection, only multipo- 
lar fields would then be possible. However, when decreasing the 
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Ekman number towards more realistic values, the simulations by 
Christensen & Aubert (2006) suggest that the dipolar window 
may persist at higher supercriticalities. In addition, since these 
stars have huge density contrasts, our investigation would then 
generally predict multipolar fields. A further exploration of the 
parameter space seems here required to clarify this point. Stanley 
& Glatzmaier (2010), for example, suggest that lower Prandtl 
number may help in creating stronger dipole fields. Considering 
radial-dependent properties (e.g. viscosity, thermal diflfusivity 
and electrical diflTusivity) is also known to have a strong impact 
on the location of the convective columns that could possibly 
help to avoid the concentration of helicity close to the equator. 
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Table 2. Results table. 



Model 




BC 




Rayleigh 


Ra/Ra^ 


/dip 




Rm 


Ro 




Ro^ 




A 






Id 


0.2 


mixed 


0.00 


8.00 X 


10^ 


9.2 


0.72 




122 


6.57 X 10" 


-3 


1.95 X 10" 


-2 


1.3 


7.13 X 10" 


-2 


2d 


0.2 


mixed 


0.00 


1.00 X 


10^ 


11.5 


0.71 




143 


7.77 X 10" 


-3 


2.55 X 10" 


-2 


2.2 


1.14 X 10" 


-1 


2m 


0.2 


mixed 


0.00 


1.00 X 


10^ 


11.5 


0.05 




180 


9.55 X 10" 


-3 


2.62 X 10 


-2 


0.6 


3.33 X 10" 


-2 


3d 


0.2 


mixed 


0.00 


2.00 X 


10^ 


23.0 


0.66 




252 


1.36 X 10" 


-2 


5.81 X 10" 


-2 


5.6 


2.24 X 10 


-1 


3m 


0.2 


mixed 


0.00 


2.00 X 


10^ 


23.0 


1.62 X 10' 


-2 


318 


1.70 X 10' 


-2 


6.26 X 10- 


-2 


2.6 


9.73 X 10- 


-2 


4d 


0.2 


mixed 


0.00 


3.00 X 


10^ 


34.5 


1.14 X 10' 


-2 


424 


r\ r\ ^ 1 /-v_ 

2.26 X 10 


-2 


8.38 X 10 


-2 


4.5 


1.42 X 10 


-1 


5d 


0.2 


mixed 


0.00 


4.00 X 


10^ 


45.9 


1.05 X 10' 


-2 


505 


2.71 X 10' 


-2 


1.05 X 10- 


-1 


6.4 


1.78 X 10- 


-1 


6d 


0.2 


mixed 


0.50 


7.00 X 


106 


5.1 


0.68 




104 


5.26 X 10" 


-3 


1.91 X 10" 


-2 


1.6 


8.68 X 10" 


-2 


7d 


0.2 


mixed 


0.50 


1.00 X 


10^ 


7.3 


0.71 




146 


7.39 X 10- 


-3 


3.28 X 10- 


-2 


4.0 


1.84 X 10- 


-1 


7m 


0.2 


mixed 


0.50 


1.00 X 


10^ 


7.3 


0.02 




185 


9.35 X 10" 


-3 


3.75 X 10- 


-2 


1.0 


4.93 X 10 


-2 


8d 


0.2 


mixed 


0.50 


1.50 X 


10^ 


10.9 


0.67 




220 


1.10 X 10" 


-2 


5.52 X 10" 


-2 


7.3 


2.83 X 10" 


-1 


8m 


0.2 


mixed 


0.50 


1.50 X 


10^ 


10.9 


0.03 




285 


1.44 X 10 


-2 


6.24 X 10" 


-2 


2.5 


r\ At . . t r\- 

9.41 X 10 


-2 


9d 


0.2 


mixed 


0.50 


1.90 X 


10^ 


13.8 


0.02 




351 


1.77 X 10" 


-2 


7.77 X 10" 


-2 


3.7 


1.18 X 10" 


-1 


lOd 


0.2 


mixed 


0.50 


2.00 X 


10^ 


14.5 


1.82 X 10" 


-2 


364 


1.83 X 10" 


-2 


8.16 X 10" 


-2 


4.0 


1.26 X 10" 


-1 


lld 


0.2 


mixed 


0.50 


3.00 X 


10^ 


21.8 


1.05 X 10' 


-2 


491 


2.47 X 10 


-2 


1.10 X 10- 


-1 


7.0 


1.78 X 10- 


-1 


12d 


0.2 


mixed 


1.00 


6.00 X 


10^ 


3.1 


0.86 




71 


3.59 X 10" 


-3 


1.55 X 10" 


-2 


0.8 


4.76 X 10" 


-2 


13d 


0.2 


mixed 


1.00 


8.00 X 


10^ 


4.1 


0.67 




117 


5.87 X 10" 


-3 


3.03 X 10" 


-2 


3.4 


1.62 X 10" 


-1 


13m 


0.2 


mixed 


1.00 


8.00 X 


10^ 


4.1 


0.70 




115 


5.81 X 10" 


-3 


2.96 X 10" 


-2 


3.2 


1.43 X 10" 


-1 


14d 


0.2 


mixed 


1.00 


1.00 X 


10^ 


5.2 


0.68 




152 


7.65 X 10" 


-3 


4.18 X 10 


-2 


5.7 


2.37 X 10" 


-1 


14m 


0.2 


mixed 


1.00 


1.00 X 


10^ 


5.2 


0.69 




149 


1A9 X 10 


-3 


A f\ A ^ ^ 'X f\- 

4.04 X 10 


-2 


6.2 


2.39 X 10 


-1 


15d 


0.2 


mixed 


1.00 


1.50 X 


10^ 


7.8 


0.52 




256 


1.29 X 10" 


-2 


1A4 X 10 


-2 


9.0 


2.83 X 10" 


-1 


15m 


0.2 


mixed 


1.00 


1.50 X 


10^ 


7.8 


0.03 




299 


1.50 X 10" 


-2 


7.51 X 10" 


-2 


3.8 


1.13 X 10" 


-1 


16d 


0.2 


mixed 


1.00 


2.00 X 


10^ 


10.3 


1.40 X 10 


-2 


390 


1.96 X 10" 


-2 


1.01 X 10- 


-1 


6.3 


1.61 X 10- 


-1 


17d 


0.2 


mixed 


1.50 


7.00 X 


10^ 


2.7 


0.81 




85 


4.29 X 10" 


-3 


2.61 X 10" 


-2 


2.0 


9.12 X 10" 


-2 


17m 


0.2 


mixed 


1.50 


7.00 X 


10^ 


2.7 


0.15 




88 


4.56 X 10" 


-3 


2.05 X 10" 


-2 


0.5 


2.26 X 10" 


-2 


18d 


0.2 


mixed 


1.50 


8.00 X 


10^ 


3.1 


0.75 




112 


5.63 X 10" 


-3 


3.57 X 10" 


-2 


3.5 


1.53 X 10" 


-1 


18m 


0.2 


mixed 


1.50 


8.00 X 


10^ 


3.1 


0.13 




126 


6.36 X 10" 


-3 


3.59 X 10" 


-2 


1.4 


5.13 X 10" 


-2 


19d 


0.2 


mixed 


1.50 


1.00 X 


10^ 


3.8 


0.63 




165 


8.29 X 10 


-3 


5.28 X 10" 


-2 


5.8 


2.03 X 10" 


-1 


19m 


0.2 


mixed 


1.50 


1.00 X 


10^ 


3.8 


0.05 




184 


9.27 X 10 


-3 


5.32 X 10" 


-2 


2.4 


7.63 X 10" 


-2 


20d 


0.2 


mixed 


1.50 


1.30 X 


10^ 


5.0 


0.05 




256 


1.29 X 10" 


-2 


7 05 X 10- 


-2 


3.7 


9.77 X 10- 


-2 


21d 


0.2 


mixed 


1.50 


1.40 X 


10^ 


5.4 


1.72 X 10" 


-2 


278 


1.40 X 10" 


-2 


7.50 X 10" 


-2 


3.9 


9.63 X 10" 


-2 


22d 


0.2 


mixed 


1.50 


1.50 X 


10^ 


5.7 


1.57 X 10" 


-2 


299 


1.50 X 10" 


-2 


8.18 X 10" 


-2 


4.5 


1.07 X 10" 


-1 


23d 


0.2 


mixed 


1.50 


2.00 X 


1 f\l 

10' 


7.6 


1.82 X 10" 


-2 


410 


2.06 X 10" 


-2 


1.17 X 10- 


-1 


8.5 


1.76 X 10- 


-1 


24d 


0.2 


mixed 


1.70 


9.00 X 


10^ 


3.1 


0.73 




134 


6.73 X 10" 


-3 


4.51 X 10" 


-2 


4.4 


1.69 X 10" 


-1 


24m 


0.2 


mixed 


1.70 


9.00 X 


10^ 


3.1 


0.07 




153 


7.71 X 10" 


-3 


4.62 X 10" 


-2 


2.0 


7 00 X 10" 

/ .V/V/ /\ XV/ 


-2 


25d 


0.2 


mixed 


1.70 


1.20 X 


10^ 


4.1 


0.06 




230 


1.16 X 10- 


-2 


6.58 X 10- 


-2 


3.4 


8.84 X 10- 


-2 


26d 


0.2 


miYPrl 


2.00 


7.50 X 


10^ 


2.2 


1.97 X 10" 


-2 


87 


4.42 X 10" 


-3 


2.31 X 10- 


-2 


0.5 


2.27 X 10" 


-2 


27d 


0.2 


mixed 


2.00 


8.70 X 


10^ 


2.5 


1.22 X 10" 


-2 


128 


6.48 X 10- 


-3 


3.80 X 10- 


-2 


1.4 


4.76 X 10" 


-2 


28d 


0.2 


mixed 


2.00 


1.00 X 


10' 


2.9 


1.72 X 10" 


-2 


165 


8.34 X 10" 


-3 


4.91 X 10" 


-2 


2.1 


6.48 X 10" 


-2 


29d 


0.2 


mixed 


2.00 


1.50 X 


\0' 


4.3 


1.37 X 10" 


-2 


294 


1.48 X 10" 


-2 


8.89 X 10" 


-2 


5.3 


1.16 X 10" 


-1 


30d 


0.2 


mixed 


2.00 


2.00 X 


10' 


5.8 


1.82 X 10- 


-2 


415 


2.08 X 10- 


-2 


1.37 X 10- 


-1 


12.8 


2.01 X 10- 


-1 


31d 


0.2 


mixed 


2.50 


1.00 X 


10^ 


2.2 


4.74 X 10 


-3 


150 


7.61 X 10" 


-3 


4.54 X 10 


-2 


2.1 


5.63 X 10" 


-2 


32d 


0.2 


mixed 


2.50 


1.50 X 


10' 


3.3 


8.63 X 10" 


-3 


282 


1.42 X 10" 


-2 


9.26 X 10" 


-2 


6.3 


1.21 X 10" 


-1 


33d 


0.2 


mixed 


3.00 


1.05 X 


10' 


2.3 


2.23 X 10- 


-3 


144 


7.28 X 10- 


-3 


3.75 X 10- 




1.3 


3.21 X 10- 




34d 


0.2 




3.00 


l.lOx 


10' 


2.4 




-3 


157 


7.92 X 10" 


-3 


4.87 X 10" 


-2 


2.4 


S 66 y 1 0" 


-2 


35d 


0.2 


mi YF»rl 

llllA.t/VJ. 


3.00 


1.20 X 


10' 


2.6 


5.60 X 10" 


-3 


184 


9.30 X 10" 


-3 




-2 


4.0 


8.27 X 10" 


-2 


36d 


0.2 




3.00 


1.30 X 


10' 


2.8 


9.43 X 10" 


-3 


211 


1 06 X 1 0" 


-2 


7.61 X 10" 


-2 


5.8 


1 OS X 10" 


-1 


37d 


0.2 




3.00 


1.50 X 


10' 


3.2 


1.84 X 10" 


-2 


261 


1.31 X 10" 


-2 


9.94 X 10" 


-2 


8.5 


1.46 X 10" 


-1 


38d 


0.2 




3.00 


2.00 X 


10' 


4.3 


8.90 X 10- 


-3 


372 


1 86 X 10- 


-2 


1.46 X 10- 


-1 


15.5 


1 68 X 10- 

X .V/O x\ XV/ 


-1 


39d 


0.2 




3.00 


4.00 X 


10' 


8.6 


1.34 X 10- 


-2 


694 


3.47 X 10- 


-2 


2 61 X 10' 

^*\J 1 r\ XV/ 


-1 


47.3 


3 18 X 10- 

mJ • I. KJ X\ X V/ 


-1 


40d 


0.6 


Stress-free 


0.01 


6.00 X 


lO 


3.5 


0.91 




45 


2.27 X 10- 


-3 


1.43 X 10- 


-2 


4.0 


1.21 




40m 


0.6 


stress-free 


0.01 


6.00 X 


10^ 


3.5 


1.06 X 10-4 


64 


3.28 X 10- 


-3 


9.78 X 10- 


-3 


0.0 


2.86 X 10- 


-3 


41d 


0.6 


stress-free 


0.01 


9.00 X 


10^ 


5.2 


0.85 




71 


3.56 X 10" 


-3 


3.15 X 10" 


-2 


0.3 


8.39 X 10" 


-2 


41m 


0.6 


stress-free 


0.01 


9.00 X 


10^ 


5.2 


2.63 X 10" 


-3 


74 


3.74 X 10" 


-3 


2.81 X 10" 




0.2 


2.72 X 10" 




42d 


0.6 


stress-free 


0.01 


1.00 X 


10^ 


5.8 


0.81 




80 


4.01 X 10" 


-3 


3.68 X 10" 


-2 


0.4 


9.89 X 10" 


-2 


42m 


0.6 


stress-free 


0.01 


1.00 X 


10' 


5.8 


2.16 X 10" 


-3 


84 


4.24 X 10" 


-3 


3.34 X 10" 


-2 


0.2 


3.09 X 10" 


-2 


43d 


0.6 


stress-free 


0.01 


2.00 X 


10' 


11.5 


0.75 




157 


7.86 X 10- 


-3 


7.02 X 10- 


-2 


3.8 


4.82 X 10- 


-1 


43m 


0.6 


stress-free 


0.01 


2.00 X 


10' 


11.5 


7.46 X 10- 


-3 


203 


1.02 X 10- 


-2 


8.98 X 10- 


-2 


1.0 


1.17 X 10- 


-1 


44d 


0.6 


stress-free 


0.01 


3.00 X 


10' 


17.3 


0.64 




241 


1.21 X 10- 


-2 


1.13X 10- 


-1 


8.5 


7.95 X 10- 


-1 


44m 


0.6 


stress-free 


0.01 


3.00 X 


10' 


17.3 


4.32 X 10- 


-3 


317 


1.58 X 10- 


-2 


1.39 X 10- 


-1 


3.0 


2.56 X 10- 


-1 
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T. Gastine et al.: Dipolar versus multipolar dynamos: the influence of the background density stratification 
Table 2. Continued. 



Model 




BC 




Rayleigh 


Ra/Ra^ 


/dip 




Rm 


Ro 




Ro^ 




A 




45d 


0.6 


stress-free 


0.01 


4.00 X 10^ 


23.0 


'^A 1 A- 

2.20 X 10 


-3 


403 


^ A'^ 1 A- 

2.02 X 10 


-2 


1.76 X 10" 


-1 


5.1 


3.67 X 10"^ 


46d 


0.6 


stress-free 


0.50 


8.00 X 


10^ 


2.6 


0.82 




58 


2.96 X 10- 


-3 


2.60 X 10" 


-2 


4.7 


9.24 X 10-1 


46m 


0.6 


stress-free 


0.50 


O A A v» 

8.00 X 


10 


2.6 


7.09 X 10- 


-3 


50 


2.53 X 10- 


-3 


2.60 X 10" 


-2 


0.1 


2.31 X 10-2 


47d 


0.6 


stress-free 


0.50 


1.00 X 


1 a6 

10 


3.2 


0.80 




73 


3.71 X 10- 


-3 


3.20 X 10- 


-2 


5.4 


8.90 X 10-1 


47m 


0.6 


stress-free 


0.50 


1.00 X 


1 n6 
10 


3.2 


2.59 X 10- 


-3 


69 


3.47 X 10- 


-3 


3.79 X 10- 


-2 


0.2 


3.57 X 10-2 


AO A 

48a 


0.6 


stress-free 


0.50 


2.00 X 


1 n6 
10 


£i A 

6.4 


0.64 




1 /C /I 

164 


8.23 X 10" 


-3 


T AA 1 A- 

/.99 X 10 


-2 


8.3 


n OA 1 A— 1 

/.80 X 10 


48m 


O.o 


stress-free 


0.50 


2.00 X 


10 


/C A 

6.4 


1.03 X 10- 


-2 


1 AC 

195 


9.77 X 10" 


-3 


1.01 X 10" 


-1 


1.2 


1 O C 1 A— 1 

1.25 X 10 


49a 


0.6 


stress-free 


A CA 

0.50 


O AA 

3.00 X 


1 a6 
10 


A C 

9.6 


A CO 

0.53 




276 


1 OA 1 A- 

1.39 X 10 


-2 


1.43 X 10" 


-1 


o o 

8.8 


C CO 1 A— 1 

6.52 X 10 


49m 


0.6 


stress-free 


A CA 

0.50 


O AA 

3.00 X 


10 


A C 

9.6 


1.10 X 10- 


-3 


O 1 o 

313 


1 CO 1 A- 

1.58 X 10 


-2 


1 CO 1 A- 

1.53 X 10 


-1 


o n 

5.1 


o on 1 A— 1 
2.8 / X 10 


50a 


0.6 


stress-free 


A CA 

0.50 


A f\f\ 

4.00 X 


10 


1 o 

lz.8 


1 OC 1 A- 

1.86 X 10 


-3 


A t A 

414 


^ AO 1 A- 

2.08 X 10 


-2 


^ AA 1 A- 

2.00 X 10 


-1 


c c 

6.5 


A AC 1 A— 1 

4.05 X 10 


51a 


0.6 


stress-free 


A CA 

0.50 


C AA 

5.00 X 


10 


1 A 

16.0 


A OC 1 A- 

9.86 X 10 


-4 


CAT 

507 


^ CC 1 A- 

2.55 X 10 


-2 


^ /I C 1 A- 

2.45 X 10 


-1 


A C 

9.5 


C AC 1 A— 1 

5.05 X 10 


52d 


0.6 


stress-free 


1.00 


1.00 X 


10^ 


1.9 


0.73 




64 


3.23 X 10- 


-3 


3.68 X 10" 


-2 


3.9 


5.72 X 10-1 


52m 


0.6 


stress-free 


1.00 


1.00 X 


10 


1.9 


1.50 X 10- 


-2 


50 


2.51 X 10- 


-3 


3.49 X 10" 


-2 


0.2 


2.21 X 10-2 


53d 


0.6 


stress-free 


1.00 


1.50 X 


10 


2.9 


0.68 




113 


5.70 X 10- 


-3 


6.98 X 10- 


-2 


3.6 


4.54 X 10-1 


53m 


0.6 


stress-free 


1.00 


1.50 X 


1 a6 
10 


2.9 


1.71 X 10- 


-2 


109 


5.47 X 10- 


-3 


6.81 X 10" 


-2 


0.6 


6.88 X 10-2 


C A A 

54a 


0.6 


stress-free 


1.00 


2.00 X 


1 a6 
10 


O A 

3.9 


0.62 




149 


7.50 X 10" 


-3 


O AO 1 A- 

8.02 X 10 


-2 


9.4 


O OA 1 A — 1 

8.20 X 10 


54m 


0.6 


stress-free 


1 AA 

1.00 


O A A 
Z.OO X 


10 


O A 

3.9 


1.55 X 10- 


-2 


173 


8.67 X 10" 


-3 


1.02 X 10" 


-1 


1 A 

1.4 


1.29 X 10-1 


CCA 

55a 


A C 

0.6 


stress-free 


1 AA 

1.00 


O AA 

3.00 X 


1 a6 
10 


c o 
5.8 


A AC\ 

0.49 




T7C 

275 


1 O O 1 A- 

1.38 X 10 


-2 


1.59 X 10" 


-1 


O A 

8.4 


c 1 A— 1 
5.65 X 10 


55m 


A C 

0.6 


stress-free 


1 AA 

1.00 


O AA 

3.00 X 


10 


c o 
5.8 


O /I C 1 A- 

3.45 X 10 


-3 


287 


1 A A ■^y t A- 

1.44 X 10 


-2 


1.59 X 10" 


-1 


A A 

4.0 


O '7 1 1 A— 1 

2. /I X 10 


en A 

5oa 


A C 

0.6 


stress-free 


1 AA 
1.00 


A f\f\ 

4.00 X 


10 


7.7 


A AA 1 A- 

4.90 X 10 


-3 


378 


1 OA 1 A- 

1.89 X 10 


-2 


O AO 1 A- 

2.03 X 10 


-1 


C A 

6.9 


O A/1 1 A— 1 

3.94 X 10 


cnA 

57a 


A C 

0.6 


stress-free 


1 AA 
1.00 


5.00 X 


10^ 


A n 

9.7 


1 0'7 1 A- 

1.37 X 10 


-3 


/I /;c 

465 


oo 1 A- 

2.33 X 10 


-2 


O /I '7 1 A- 

2.4/ X 10 


-1 


1 A 1 

10.1 


C AO 1 A— 1 

5.03 X 10 


58d 


0.6 


stress-free 


1.50 


1.20 X 


10 


1.5 


0.66 




58 


2.93 X 10- 


-3 


3.56 X 10" 


-2 


2.9 


3.79 X 10-1 


58m 


0.6 


stress-free 


1.50 


1.20 X 


10 


1.5 


0.06 




46 


2.32 X 10- 


-3 


3.03 X 10" 


-2 


0.2 


3.03 X 10-2 


59d 


0.6 


stress-free 


1.50 


1.50 X 


10 


1.8 


0.66 




89 


4.47 X 10" 


-3 


6.08 X 10- 


-2 


3.5 


3.77 X 10-1 


59m 


A 

0.6 


stress-free 


1 CA 

1.50 


1.50 X 


1 a6 
10 


1.8 


A A'7 

0.07 




81 


A AA 1 A- 

4.09 X 10 


-3 


5.64 X 10" 


-2 


A C 

0.5 


A on 1 A— 2 

4.8 / X 10 


oOd 


A 

0.6 


stress-free 


1 CA 

1.50 


2.00 X 


1 n6 
10 


o c 

2.5 


A AC\ 

0.49 




1 /I 

146 


7.33 X 10" 


-3 


1.02 X 10" 


-1 


3.2 


2.54 X 10-1 


oOm 


A /C 

0.6 


stress-free 


1 CA 

1.50 


2.00 X 


10^ 


o c 

2.5 


A A/1 

0.04 




1 /I o 

14z 


7.12 X 10" 


-3 


A O O 1 A- 

9.28 X 10 


-2 


1 o 

1.3 


1.05 X 10-1 


nt A 

old 


A C 

0.6 


stress-free 


1 CA 

1.50 


2.50 X 


10^ 


o 1 
3.1 


1 /I C 1 A- 

1.45 X 10 


-2 


zOz 


1 A1 1 A- 
1.01 X 10 


-2 


1.25 X 10" 


-1 


o o 

3.2 


O O'T 1 A— 1 

2.2/ X 10 


n^A 

ozd 


A C 

0.6 


stress-free 


1 CA 

1.50 


3.00 X 


10^ 


3.7 
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